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ABSTRACT 

RW Aurigae (RW Aur) is a binary star system with a long molecular arm trailing 
the primary star. Cabrit et al. (2006) noted the resemblance between this extended 
structure and the tidal arm stripped from the primary star in the simulations of star- 
disc encounters by Clarke & Pringle (1993). In this paper we use new hydrodynamical 
models and synthetic observations to fit many of the parameters of RW Aur. Using 
hydrodynamic models we find that the morphological appearance of RW Aur can be 
indeed explained by a tidal encounter with the secondary star. We reproduce all the 
major morphological and kinematic features of the system. Using radiative transfer 
calculations, we find that synthetic CO and dust continuum observations of our hydro- 
dynamic models agree well with observations. We reproduce all the main features of 
the line profiles, from emission fluxes to the optical depth of the different components 
of the system. The agreement between observations and simulations thus lends strong 
support to the hypothesis of a tidal encounter scenario. Finally, we propose a possible 
solution for the origin of the dimming of the primary star observed in 2010/2011 by 
Rodriguez et al. (2013). 

Key words: accretion, accretion discs - protoplanetary discs - binaries:general - 
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1 INTRODUCTION 

RW Aurigae is a young stellar system comprising two 
roughly solar mass stars at a projected separation of 1.5'/ 
i.e. 200 AU at a distance d ~ 140 pc (cfr. Hipparcos cat¬ 
alogue, van Leeuwen 2007). The primary is a Classical T 
Tauri Star which is characterized by an unusually high ratio 
of accretion rate (2 — 10 x lO^^M© yr“^, Basri & Bertout 
1989; Hartigan, Edwards & Ghandour 1995) to disc mass 
(~ 1 O~^M 0 as derived from mm continuum dust emission, 
Osterloh & Beckwith 1995). 

The most notable feature of the RW Aurigae system is 
the 600 AU long arm structure trailing from the disc of RW 
Aur A which was imaged by Cabrit et al. (2006) in^^CO(J = 

2 — 1 ) and ^^CO( J =1 — 0) with the IRAM(Institut de Ra- 
dioastronomie Millimtrique) Plateau de Bure Interferome¬ 
ter. Cabrit et al. (2006) noted the resemblance between this 
structure and the tidal arm stripped from the primary star in 
the simulations of star-disc encounters by Clarke Sz Pringle 
(1993). Although these simulations were not designed to fit a 
particular observation, the superficial resemblance between 
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the morphology of RW Aur and the prograde encounter of 
Clarke & Pringle (1993) was striking. It is however unclear 
without detailed modelling whether simulations can fit all 
the parameters of the system, particularly the gas kinemat¬ 
ics and the available radial velocity and proper motion data 
for the two stars. Moreover Bisikalo et al. (2012) presented 
evidence from proper motion data that the encounter might 
in fact be retrograde in which case the Clarke & Pringle 
(1993) models would predict a morphology unlike that ob¬ 
served in RW Aur. Bisikalo et al. ( 2012 ) have on this basis 
proposed an alternative model involving circumbinary disc 
accretion. In the detailed modelling that we present below 
we shall attempt to fit all the available kinematic constraints 
through a star-disc encounter (i.e we assume that the system 
is either highly eccentric or marginally unbound. Eor justi¬ 
fication of this assumption, see Section 3.3), but note that, 
since the errors on the proper motion data in Bisikalo et al. 
( 2012 ) are large, we shall place more emphasis in fitting the 
radial component of the stars’ relative space motions on the 
plane of the sky. 

Interest in the RW Aurigae system has recently been 
revived by the photometric observations of Rodriguez et al. 
(2013). Their optical monitoring detected a 2 magn dim¬ 
ming event lasting 180 d which the authors intriguingly in- 
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terpreted in terms of occultation of RW Aur A by a por¬ 
tion of the tidal arm. We shall investigate whether our best 
model for the encounter involves a geometry that would per¬ 
mit such an event. 

In this paper we use new hydrodynamical models and 
synthetic observations to fit many of the parameters of RW 
Aur. We aim to understand the geometry of the binary 
system, the mechanism/geometry responsible for its mor¬ 
phological appearance, other observed properties (such as 
molecular line profiles) and the recent dimming event found 
by Rodriguez et al. (2013). The rest of this paper proceeds 
as follows. In Section 2 we describe the range of observa¬ 
tional constraints of the RW Aur system which we shall use 
to constrain our models. In Section 3 we outline the codes 
used in our investigation (the PHANTOM^ smooth particle 
hydrodynamics code which is used for the hydrodynamical 
simulations and the TORUS^ Monte Carlo radiation transport 
code which is used to produce synthetic observations from 
the results of the hydrodynamic models) and also the initial 
conditions and parameter space surveyed. In Section 4 we 
present and discuss our best-fitting model solution and the 
degree to which various system parameters are constrained 
by the observations. In Section 5 we present the detailed 
line profiles and continuum emission map for our best-fitting 
model. Section 6 discusses a possible origin of the ‘dimming 
event’ in the light curve of RW Aur, and whether it could 
be associated with occultation by the tidal arm. Section 7 
summarizes our conclusions. 

2 OBSERVATIONS OF RW AUR 

2.1 The disc around star A 

The ^^CO( J = 2 — 1) emission from close proximity to RW 
Aur A shows a perfect point symmetry between its red- 
shifted and blueshifted lobes. Furthermore, the radial ve¬ 
locity gradient is perpendicular to the optical jet. Both of 
these features suggest an axisymmetric geometry, i.e. RW 
Aur A most likely has a rotating disc (Cabrit et al. 2006). 
The sky inclination of the disc was constrained using the ra¬ 
tio of the radial and tangential velocities of emission knots 
in its jet (Lopez-Martm, Cabrit & Dougados 2003). Alter¬ 
natively, the inclination was also constrained by the radial 
velocity of the jet itself. Combining these, the sky inclina¬ 
tion was estimated to be in the range of 45 — 60° (Woitas 
et al. 2002). 

The observed flux densities in ^^CO(J = 2— 1) and 
^^CO( J =1 — 0) has a ratio of about 4 (see fig. 4 by Cabrit 
et al. 2006). The ratio of 4 is the square of the frequency 
ratio between these two lines. This suggests that the disc 
is optically thick in which case the Rayleigh-Jeans limit 
predicts that Bu — 2kTv^!(?. Beckwith h Sargent (1993, 
hereafter BS93) predicted that for an optically thick, Keple- 
rian, vertically isothermal disc of radial temperature profile 
T oc the molecular line profile scales as towards 

the wings of the line profile. The disc line profile observed 
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by Cabrit et al. (2006) suggests that q ~ 0.73. The same pa¬ 
rameter q was also independently estimated as ~ 0.57 from 
the infrared spectral energy distribution (Osterloh & Beck¬ 
with 1995). The BS93 model also predicted that the peak 
flux density of disc line profiles occurs at the velocity of the 
disc’s outer radius: Vpeak = sin (z)(GM/Rout)°'^- The ob¬ 
served peak for A’s disc (~ 3.5 km/s) then constrains Rout 
to be in the range of 50 — 77 AU respectively for an inclina¬ 
tion of 45 — 60°. Cabrit et al. (2006) noted that stellar mass 
and outer disc radius are degenerated when considering the 
observed peak velocity. Therefore they argued that the bet¬ 
ter approach is to compare the shape (rather than absolute 
peak velocity) of observed and theoretical line profiles in 
BS93; they claimed Rout to be 40 — 57 AU. 

Cabrit et al. (2006) estimated the temperature of the 
disc at its outer radius with two methods: 

1) Measuring the peak flux density of the disc line profile 

and using the BS93 model which predicted that the peak 
flux density scales as F(X) ~ d~‘^; 

2) Using the brightness temperature since the disc is 
optically thick and in the Rayleigh-Jeans regime. 

The two methods consistently gave a temperature of 60 — 
107 K (the range corresponding to different sky inclination) 
at the outer radius of the disc. 

Cabrit et al. (2006) further estimated the column den¬ 
sity of the disc. The fact that the disc is optically thick 
in ^^CO and undetected in ^^CO placed lower and upper 
limits on the column density respectively. These CO col¬ 
umn density limits were converted to total column densites 
using standard ^^CO and ^^CO abundances in the local in¬ 
terstellar medium. Assuming that the disc surface density 
scales as 1/R, Cabrit et al. (2006) derived an estimate of 
the total mass of the disc of 10“^ — lO^^M©. This is consis¬ 
tent with the mass estimate by Osterloh & Beckwith (1995) 
(3 X lO^^M©), but not with the higher estimate by Andrews 
& Williams (2005) (4 x lO^^M©). 

2.2 The tidal arm 

^^CO molecular emission from the tidal arm is entirely red- 
shifted (with the peak velocity at 3.2 —3.8 km/s with respect 
to star A, Cabrit et al. 2006). The observed radial veloc¬ 
ity is larger than the escape velocity at the corresponding 
distance, indicating that the motion of the tidal arm is un¬ 
bound. The tidal arm is also optically thick as the observed 
^^CO( J = 2 — 1)/ ^^CO(J =1 — 0) flux density ratio is ~ 4. 
The tidal arm extends about 4 (~ 600) AU at the distance 
of 140 pc. 


2.3 The molecular complex around star B 

^^CO molecular emission around star B has a broad and 
asymmetric peak blueshifted by 3 — 6 km/s with respect 
to star A. The asymmetric peak cannot be attributed to 
an undisturbed rotating disc. Given that the tidal interac¬ 
tion hypothesis is favoured, the molecular complex around 
B most likely consists of material captured during the tidal 
encounter. The ^^CO(J = 2 — 1)/ ^^CO(J =1 — 0) ra- 
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Figure 1. A schematic illustrating the polar and azimuthal an¬ 
gle used to define the mutual inclination between the angular 
momenta of the disc and the orbit of the secondary. 

tio was observed to be ~ 5.7 suggesting that the molecu¬ 
lar complex around B is marginally optically thick (Cabrit 
et al. 2006). Note, however, that the signal to noise of the 
^^CO(J = 1 — 0) emission line is very low, we therefore do 
not use this line as a constraint to discuss our model. 


2.4 Proper motion of the system 

Bisikalo et al. (2012) compiled the position angles of star 
B and the angular separations between star A and B over 
the past 70 yr. Linear regression analysis showed that the 
projected AB separation was increasing at 0.002 / yr. 70 yr 
is short compared to the Keplerian orbital period (~ 2000 
yr) at the observed separation (~ 1.5 , or 200 AU for 140 
pc away), thus justifying a linear approximation for the ve¬ 
locity. They similarly claimed that the position angle of star 
B was increasing at a mean rate of 0.02°/yr, which would 
imply that the orbital motion of the pair was anti-clockwise 
as projected on to the plane of the sky. This conclusion is 
however not firm given that the claimed increase is based on 
a fit involving data from the 1940s whose radial component 
was discrepant with the recently observed trend of increas¬ 
ing separation. Moreover this fit to the tangential data is 
in the opposite sense to the apparent trend observed dur¬ 
ing the 1990s (although the uncertainties during this period 
were large). We therefore only constrain our model such that 
the sky-projected separation of star A and B is increasing 
at the current epoch. Future proper motion observations are 
needed to decide whether the motion of the two stars pro¬ 
jected into the plane of the sky is clockwise or anti-clockwise. 


2.5 The dimming event 

Rodriguez et al. (2013) reported a long (~ 180 d) and deep 
(2 mag) dimming of RW Aur from 2010 September until 
2011 March. Typical photometric variations of CTTS (Clas¬ 
sical T Tauri Stars) only last for days and weeks with a depth 
< 1 mag. This is usually associated with circumstellar ex¬ 
tinction or hotspots on the accretion disc that rotate into 
and out of the observer’s view. Rodriguez et al. (2013) noted 


that there was no similar dimming event^ for the system in 
the past 50 yr and proposed that the dimming may be due 
to the occultation of star A by the leading edge of the tidal 
arm. They further estimated that the occulting body was 
moving with a transverse velocity of 0.86 — 2.58 km/s (di¬ 
viding the stellar radius of star A by the ingress time of the 
dimming). The size of the occulting body was estimated to 
be ~ 0.09 —0.27 AU using the product between the duration 
of the dimming and the derived velocity. The large range in 
the estimates of transverse velocity and size are due to the 
large uncertainties on the ingress time-scale, approximately 
between 10 and 30 d. 

More recently, Antipin et al. (2015) have again detected 
RW Aurigae in the ‘low state’ (i.e. 3-4 mag less than what 
observed by HST in the mid-90’s. White & Ghez 2001). 


2.6 Summary of observational constraints 

Among the various observations, we identified the following 
key criteria to gauge how closely our simulations describe 
the observations. These features were more reliably con¬ 
strained by observations; they are also easily quantifiable 
in simulations: 

1) The size of the system: Cabrit et al. (2006) reported 
that the tidal arm extended about 4 (~ 600 AU at a 
distance of 140pc); the projected AB separation was ~ 1.5 
(~ 200 AU). We aim to reproduce the size of the system 
as observed. However we argue that it is more important 
to fit the relative ratio between the dimensions of various 
structures of the system. This is because the relative ratio 
is independent of the error in the distance measurement (a 
parallax of 7.08 ± 0.71 mas corresponding to a distance of 
140 ± 14 AU Basri k Bertout 1989). 

2) The relative orientation of the A disc and star B. We 
gauged the relative orientation using position angles (PA). 
For the disc, its position angle is defined as the angle 
between its radial velocity gradient and north. For star B, 
PA is defined as the angle between the vector joining star 
A and B and North. According to Cabrit et al. (2006), the 
PA of the disc is ~ 40°; the PA for B is ~ 105° (or 255° if 
using the convention of Bisikalo et al. 2012). 

3) Sky inclination of the disc: the sky inclination of the disc 
was estimated to be in the range of 45 — 60° (Woitas et al. 
2002; Lopez-Martm, Cabrit k Dougados 2003). 

4) Proper motion: as reported by Bisikalo et al. (2012) 
the angular separation between star A and star B was 
increasing at a rate of 0.002''/ yr. 

5) The line profiles: the radial velocity of the peak, the peak 
flux density and the optical depth. 


^ Note that the long term light curve does show a number indi¬ 
vidual photometric data points at the ‘low’ level recorded in the 
2011 event but a well resolved eclipse light curve has not been 
detected hitherto. 
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Figure 2. The gas morphology of the best-fitting model in the 
frame of reference comoving with star A. The x-y plane indicates 
the plane of the sky. Crosses indicate SPH particles, which were 
selected randomly every other 10 particles. The black stars show 
the position of Star A and B. Colour coding is linear in radial 
velocity (km/s). Note, for direct comparison with Fig 3(c), the 
heliocentric velocity (15.87T 0.55 km/s, Cabrit et al. 2006) of the 
star A should be added. 

3 NUMERICAL METHOD 

3.1 Hydrodynamics models 

We model the stellar fly-by using the smooth particle hy¬ 
drodynamics (SPH) code PHANTOM (see e.g. Lodato & Price 
2010; Price & Federrath 2010; Price 2012; Facchini, Lodato 
& Price 2013; Nixon, King & Price 2013). We assume that 
self-gravity in the disc is negligible and that the gravitational 
potential is dominated by the stars. We employ a locally 
isothermal equation of state and adopt standard artificial 
viscosity parameters {a = 0.1;/3 = 2;cr = 0.1) according to 
the presciption of Morris & Monaghan (1997). 

3.2 Synthetic observations 

We use the radiation transport and hydrodynamics code 
TORUS to compute synthetic observations in this paper. 
TORUS is predominantly a Monte Carlo radiation transport 
code used for calculating line and continuum transport (e.g. 
Harries 2000; Kurosawa, Harries & Symington 2006) but 
also features photoionization (Haworth, Harries & Acreman 
2012), non-LTE molecular line transfer (Rundle et al. 2010) 
and hydrodynamics (Haworth & Harries 2012). 

We map the result of the SPH simulation on to the 
TORUS grid using the technique described in Acreman, Har¬ 
ries & Rundle (2010). We do not explicitly derive the tem¬ 
perature using radiative transfer, but rather assume that 
the temperature distribution in the model is a function of 
cylindrical distance from the primary (see below; equation 
1), as in Cabrit et al. (2006). For molecular line transfer 
calculations we solve the equations of non-LTE statistical 
equilibrium and produce synthetic data cubes in the man¬ 
ner described in detail by Rundle et al. (2010). 

Since we are comparing with the observations from 
Cabrit et al. (2006) we consider the ^^CO(J = 2—1) and 
(J = 1 — 0) lines and assume the same parameters, i.e. the 
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Figure 3. Fig. 1 from Cabrit et al. (2006), showing the inter¬ 
ferometric data of PdBl. From their caption: “a) 1.3 mm contin¬ 
uum map of the RW Aur system. Contour levels are 3, 6, 12, 24, 
48cr (Icr = 0.47mJy/beam). Crosses indicate the peak positions, 
which closely correspond to RW Aur A and B. The large arrow 
shows the direction of the blueshifted jet. b) ^^CO(J = 2 — 1) 
emission blueshifted and redshifted with respect to the RW Aur 
A systemic velocity. Contour spacing is 3cr (0.12 Jy km s“^ ), 
truncated to 18f7 for clarity, c) and d) Heliocentric centroid ve¬ 
locity and line width (i.e. first and second order moments) of the 
i2CO(J = 2-1) line.” 


^^CO abundance relative to hydrogen is 8 x 10“^ and the 
system is at a distance of 140 pc from the observer. We 
assume a small micro turbulence value of 0.03 km/s. The 
data cubes comprise 401 x 401 pixels and 36 velocity chan¬ 
nels spanning —7.19 to 7.63 km/s. This velocity range and 
number of channels is chosen to match the IRAM Plateau 
de Bure Interferometer, the instrument used for the obser¬ 
vations in Cabrit et al. (2006) against which we compare 
our models. The beam size is chosen to match that of the 
IRAM Plateau de Bure Interferometer (2.42" x 1.55" and 
0.89" X 0.58" for the J = 1 — 0 and J = 2 — 1 transitions 
respectively). We convolve the resulting data cubes with a 
2-D Gaussian with a width corresponding to the semiminor 
axis of the two beams using ACONVOLVE from CIAO v 4.5 
(Fruscione et al. 2006). Finally, the line profiles of the syn¬ 
thetic observations are obtained by integrating over ellipses 
that match the ones used by Cabrit et al. (2006). 

We also compute a 1.3 mm continuum image using 
Monte Carlo radiative transfer, algorithmic details of which 
are given in e.g. Harries (2000) and Harries et al. (2004). We 
discuss the assumed dust distribution in Section 5.2. 
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3.3 Initial conditions 

We model the encounter geometry shown in Fig 1. in which 
a spherical polar system (with the z-axis along the normal 
to the plane of the secondary’s orbit, and the x-axis point¬ 
ing towards the apocentre) is used to define the coordinates 
{0, 0) of the disc normal (where the separation vector of the 
two stars at pericentre is located at 0 = tt). We populate 
particles in a Keplerian disc around the primary according 
to a prescribed surface density law (see below) and with 
a Gaussian distribution normal to the disc plane which is 
the hydrostatic equilibrium solution set by the imposed lo¬ 
cal temperature (see below). For reasons of computational 
economy we set an inner radius of 5.9AU and do not con¬ 
tinue to integrate particles that stray within this radius. 

We argue that the presence of a single well-defined tidal 
feature in the observations suggests that the disc of RW 
Aur A has been subject to either a one-off encounter or to 
an encounter that is well separated in time from any pre¬ 
vious interactions. Thus we consider it likely that the orbit 
is either of high eccentricity or unbound. It is unlikely that 
the encounter is highly hyperbolic since this would imply a 
relative velocity at infinity that exceeds typical stellar ve¬ 
locity dispersions in star forming regions. We thus expect 
that the orbit - whether mildly bound or mildly unbound - 
is likely to be close to parabolic and hence model this sit¬ 
uation here.This is consistent with the low relative radial 
velocity between the two stars reported by White & Hillen- 
brand (2004). The masses of RW Aur A and B are set to 
1.4Mo and 0.9Mq respectively (Ghez, White & Simon 1997; 
Woitas, Leinert & Kohler 2001). 

We set the surface density profile to S oc 1/R and adopt 
a power law temperature profile as a function of cylindrical 
distance from the primary parametrized according to 

T = To(R/Ro)“" (1) 

The radial index, q has been estimated to be in the 
range 0.5 — 0.73 by various authors {q = 0.57 and q = 0.73 
Osterloh & Beckwith 1995; Gabrit et al. 2006, respectively). 
We adopt a temperature around 60 K at R = 57 AU in line 
with the estimate of Gabrit et al. (2006) for an inclination 
of the disc i = 60°, and vary the value of q (which affects 
the pressure scale length of the outer disc and has a modest 
effect on the morphology of the tidal tails produced). 

The simulation is initiated when the stellar separation 
is > 10 pericentre radii in order to ensure that the disc is 
not violently perturbed by the sudden introduction of the 
companion in its vicinity. 


3.4 Constraining other parameters 

After setting up our model with the plausible initial condi¬ 
tions discussed in the previous section, in principle one could 
constrain the various other orbital parameters by running 
the hydrodynamic simulation and synthetic observations in 
conjunction with a parameter study tool (e.g. Markov chain 
Monte Garlo). However, such a parameter study is unfeasi¬ 
ble in our case: conventionally, 10® particles are needed for 
sufficient resolution in such a hydrodynamic simulation. In 
conjunction with the simulated observations, this would take 
several days of computation per model. Moreover there are 


at least 10 remaining parameters in the model and so the 
time required for such an approach would be prohibitive. 

Instead we adopted the following strategy. We initially 
ran our hydrodynamic code with 10^ test particles on a 
coarse grid in the 10-D parameter space. We then focused 
on the region that is most consistent with various observa¬ 
tional constraints. Parameters were then tuned until simu¬ 
lated values converged to observed values within respective 
uncertainties. Finally we repeated the same model with 10® 
test particles and appropriate scalings (e.g. viscosity scales 
with number of particles) to reveal small-scale details of the 
system. 

We remind the reader that although we refer to this fi¬ 
nal model as the best-fitting one, in reality this is only the 
best-fitting model among the few dozen models we tested. 
Given the uncertainties in the observations (Section 2) and 
the various shortcomings of our simulations, the statistical 
errors introduced from our sampling strategy are likely to 
be outweighed by other factors (see discussion in Section 5). 
We nevertheless argue that the good level of agreement be¬ 
tween simulations and observations in the models performed 
is already sufficient to demonstrate that the tidal encounter 
hypothesis is indeed consistent with the existing observa¬ 
tions. 


4 RESULTS OF HYDRODYNAMICS MODELS 

In this section, we first consider the parameters of the en¬ 
counter that affect the gross properties of the flow (i.e. the 
morphology and kinematics of the tidal tail and disc). It 
is sufficient for such an analysis to be conducted with the 
hydrodynamic code alone in the limit of optically thin emis¬ 
sion. We fix the stellar masses and eccentricity of the orbit 
for all simulations and for each simulation we adopt given 
values of 0 and 0, of pericentre distance, disc outer radius, 
surface density and temperature normalization, and temper¬ 
ature power-law index. In Section 5, we consider whether a 
plausible initial disc model can then reproduce line fluxes 
and line profiles by using synthetic observations for a range 
of output times and viewing angles (constrained so as to give 
the correct projected separation of the two stars on the sky). 
The viewing angle is parameterised in terms of the direction 
of the line-of-sight vector identified by the polar coordinates 
^obs and 0obs in the coordinate frame shown in Fig. 1. 

4.1 Sensitivity of results to parameter choice 

To determine which model is most consistent with observa¬ 
tions, we vary each parameter around the initial estimates 
obtained from the pilot simulations. We then identify the 
system observables that are most dependent on each pa¬ 
rameter in order to constrain that parameter. 

The choice of periastron separation (70 AU) was mainly 
decided by fitting observed length and width of the tidal tail 
both of which increase with pericentre separation (see Fig. 
5a). The length and width of the tidal tail is also dependent 
on time of observation. However this degeneracy can be bro¬ 
ken, since the choice of periastron separation must also give 
a consistent radial velocity of material in the tidal tail which 
is similar to the escape velocity at periastron (see Fig. 8). 
Moreover, since the disruptiveness of the encounter depends 
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Figure 4. Evolution of the tidal encounter in the plane of the sky, shown with four different snapshots. The bottom right panel shows 
the structure as seen today. From top left to right bottom the snapshots are taken at —589, —398, —207 and 0 yr from today, respectively. 
The surface density is logarithmically scaled, and was chosen such that the final snapshot resembles the observations by Cabrit et al. 
(2006). Note that a stream of material linking star A and star B is present, even though it does not appear in the bottom right panel, 
due to the choice of surface density threshold (compare this final snapshot with the right-hand panels of Fig. 11). 


Table 1. Quantitative comparison of best-fitting model outputs 
and observed values. “ From Cabrit et al. (2006) assuming d = 140 
pc where necessary. ^ Lopez-Martfn, Cabrit & Dougados (2003), 
Woitas et al. (2002). ^ Bisikalo et al. (2012). 


Observed features 

Observed values 

Best fit simulation 

values 

Length of tidal arm 

600 AU (4")“ 

650 AU 

AB separation on sky 

200 AU 

200 AU 

PA of disc (°) 

40“ 

47.2 

PA of AB (°) 

105“ 

105.9 

Inclination of disc (°) 

0 

1 

10 

64.0 

Rate of change of 

0.002"yr“^'' 

0.003"yr“^ 


on the ratio of the periastron separation to the outer disc 
radius, the observed scale of the residual disc constrains the 
initial disc radius (60 AU). 

Another important quantity that could be used to con¬ 
strain the orbital parameters of the flyby is the disc outer 
radius, after truncation. Since Clarke & Pringle (1993) other 
groups have worked on the angular momentum and energy 
transfer mechanisms occurring with a star-disc encounter 
(e.g. Kobayashi & Ida 2001; Munoz et al. 2015). Recently, 


Breslau et al. (2014) focused on the sizes of protoplane¬ 
tary discs after such encounters using numerical simula¬ 
tions, where the discs are described by test particles (hy¬ 
dro effects are neglected). They model the encounter via 
prograde coplanar parabolic orbits. They obtain the follow¬ 
ing prescription for the outer radius of the truncated disc: 
Rout = 0.28RperiTnj“2°‘^^, where mi 2 = M 2 /Mi. However, 
note that they did not explore the effects of the misalign¬ 
ment between the disc and the plane of the parabola. Recent 
results by Lubow, Martin & Nixon (2014) analysed how the 
tidal truncation radius estimates are affected by the mis¬ 
alignment between a disc and an external binary (Artymow- 
icz & Lubow 1994). Their results indicate that the external 
torque decreases quite rapidly with inclination, leading to a 
larger circumstellar disc. Similarly, we expect the same kind 
of behaviour for a flyby interaction, even though in this case 
resonances do not come into play. Finally, the analytic pre¬ 
scription by Breslau et al. (2014) strongly depends on the 
definition of the outer radius. They define it as the point 
where the surface density profile presents a maximum in its 
radial gradient. However nothing prevents an observer from 
detecting regions of the disc outside this radius, where the 
material will be less bound and in more elliptical orbits (see 
again Breslau et al. 2014). To conclude, the outer radius is 
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Figure 5. In each of the diagrams, one parameter of the model is varied at a time and plotted against the observational feature that is 
most dependent on it. Shaded areas show observed values within Icr. Comparison with observed values helps to constrain the parameter 
being varied (see Table 2 for the best-fitting model). Top left: the pericentre against the length and width of tidal tail; top right: the 
polar angle 0 between the orbital and disc angular momenta against the positional angle of the captured material on star B; centre left: 
the azimuthal angle (j) between the orbital and disc angular momenta against the positional angle of star B; centre right: the time of 
observation against the length and width of tidal tail; bottom: the power index q of radial temperature profile against the aspect ratio 
(width/length) of the tidal tail. 
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Table 2. Parameter choices that best fit the observations and the effect of parameter variation on system properties. 


Parameter varied Best fit value Effect 


Initial disc radius 
Periastron separation 
Mutual inclination of disc and binary planes 


60 AU 
70 AU 

0 = 18°, (j) = 170° 


Viewing angle 
Time of observation 


= 114° 

l.lPperi after pericentre 


Radial Temperature profile 
Initial Disc Mass 


r = 230K(^) 
0.00167Mo 


Post-encounter disc size and tidal tail morphology 
Morphology and kinematics of tidal tail 
0 determines the orientation of tidal tail and 
the captured material on star B; 

(j) determines the velocity and positional angle of star B 
Morphology in plane of sky 

Kinematic signatures and ratio of AB separation to 
extent of tidal arm 

Power law index determines aspect ratio of tidal tail 
Optical depth of different regions of the system 


not the best parameter to be used to constrain a priori the 
orbit of the flyby. We will check a posteriori that what we 
obtain is in reasonable agreement with the observations (see 
Section 5). 

The mutual inclination of the disc and the binary orbit 
(0) is the main determinant of the morphology of the sys¬ 
tem. While the azimuthal angle 0 determines the relative 
position and separation between star A and star B, the po¬ 
lar angle 0 determines the size and orientation of the tidal 
arm. The tidal arm is most prominent and most massive for 
a prograde coplanar tidal encounter. As 0 increases, the tidal 
arm decreases in size and mass; it also becomes warped with 
respect to the original disc plane. If 0 is larger than 90°, no 
extensive tidal arm is produced. 

We therefore require a prograde encounter. We can ob¬ 
tain a clue about the likely orientation of the disc of A on 
the sky from the observation of its jet which is directed to 
the SE. Since the jet must be visible on the observer’s side 
of the disc, it suggests that the NW side of the disc is nearer 
to the observer. The fact that the NE lobe of A’s disc is 
blueshifted then suggests that the rotation of the disc is 
clock-wise as seen from the earth. In the successful (rather 
close to coplanar prograde) model, we have found a viewing 
angle that well reproduces the observed blueshift of the ma¬ 
terial centred on B. However, the predicted proper motion 
in the plane of the sky is then eloek-wise in contrast to that 
claimed by Bisikalo et al. (2012). We have discussed above 
how it is desirable to better constrain the motion in the 
plane of the sky through future proper motion observations. 

The power index q of the adopted radial temperature 
profile determines the rate at which temperature falls as a 
function of radius. A smaller value of q results in a hot¬ 
ter and thus thicker outer region of the undisrupted disc. 
Subsequently, a thicker outer disc produces a thicker tidal 
tail after the tidal interaction. Therefore, the value of q is 
constrained via the observed ratio between the width and 
length of the tidal arm. 

4.2 The best fitting model 

Eig. 2 shows the gas morphology of the best-fitting model, 
the parameters of which are listed in Table 2 (compare Eig. 
2 with Eig. 3, where we portray Eig. 1 by Cabrit et al. 2006). 
The gas is rendered as individual SPH particles which are 
colour coded according to line of sight velocity. Eig. 4 shows 
the evolution of the fly-by with rendered plots. 



1.1 2 2.8 3.7 4.5 5.4 6.2 7.1 7.9 


Figure 6. The integrated ^^CO( J = 2 — 1) emission map of RW 
Aur after correcting for beam dilution. Colour coding is linear in 
Jy beam“^ km/s. 


This model successfully reproduced the following ob¬ 
served features of the system: 

• The northeast lobe of the disc is blueshifted with ra¬ 
dial velocity of ~ —3 km/s with respect to Star A while the 
southwest lobe of the disc is redshifted with radial velocity 
of r\j 3 km/s. 

• The tidal arm has the observed spiral shape with the 
correct handedness. It is entirely redshifted with a radial ve¬ 
locity of ~ 3 km/s which is larger than the escape velocity at 
the corresponding distance. The tidal arm is thus unbound, 
as observed. 

• The molecular complex on B is blueshifted with respect 
to star A. The shape of the captured material around B is 
similar to that observed with its major axis being tilted 
slightly to the east of north. This good agreement suggests 
that the gas observed around B may indeed be material 
captured from A’s disc (even if B also originally possessed a 
small scale disc of its own). 
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Comparing the observation and simulation more quanti¬ 
tatively, Table 2 lists the various quantifiable morphological 
features. Our simulations matches the observations in nearly 
all respects. The strikingly good agreement demonstrates 
that the system’s morphology is most likely associated with 
a tidal interaction. 


5 RESULTS OF SYNTHETIC OBSERVATIONS 
5.1 Line profiles 

We took snapshots from our hydrodynamic models at ~ 570 
yr (^obs ~ 1-lPperi, whcic Pperi IS thc Kcplcrian period at the 
pericentre) after pericentre and post processed them with 
TORUS to produce synthetic data cubes (c.f. section 

3.2). Fig. 8 shows line profiles from subsets of these data 
cubes after applying convolution accounting for the IRAM 
beam size, for the best-fitting model. 

The synthetic observations allow us to constrain the 
temperature normalization and the initial mass of the disc 
Md, since these two quantities regulate the flux density nor¬ 
malization and the optical depth of the emission lines. In 
the case that the line is optically thick, the ratio ^^CO(J = 
2 — 1)/ ^^CO(J =1 — 0) should be ~ 4 in Rayleigh-Jeans 
limit. We have therefore multiplied the ^^CO(J = 1 — 0) by 
4 for the regions of the tidal arm and the primary disc, to 
demonstrate whether these components of the system are 
optically thick. We do not consider the ^^CO( J =1 — 0) line 
for the gas around star B, since this line has been marginally 
detected by Cabrit et al. (2006), as discussed above. The ra¬ 
dial velocities were corrected to heliocentric frame by adding 
the mean heliocentric velocity of star A (15.87 ± 0.55 km/s, 
Cabrit et al. 2006). 

First, we consider the impact of temperature normal¬ 
ization. By using the temperature profile as estimated by 
Cabrit et al. (2006) (T(P) = 450 K x (P/AU)“°‘^), the most 
outstanding discrepancy is that the simulated flux densities 
are 3—5 times higher than the observed flux densities. We ar¬ 
gue that the most likely cause is that the limitation of BS93 
model led to Cabrit et al. (2006) overestimating the tem¬ 
perature of the disc. Although BS93 assumed local thermal 
equilibrium (LTE), this assumption may not be appropriate 
for RW Aur. TORUS is capable of treating both LTE and 
non-LTE, allowing us to test this. The comparison in Eig. 
7 demonstrates that the resultant flux densities are indeed 
about 2 — 3 times lower in the LTE case. 

This indicates that in order to explain the observed flux, 
Cabrit et al. (2006) would have overestimated the temper¬ 
ature by assuming an LTE model. Eurthermore, if we use 
a lower temperature, more consistent with theoretical ex¬ 
pectation from stellar irradiation and accretion, (T = 30K 
at 40 AU, Chiang & Goldreich 1997), the resultant theo¬ 
retical emission is more compatible with the observed inte¬ 
grated flux densities. Note that lower temperatures are also 
indicated by Infrared SED fitting by Osterloh Sz Beckwith 
(1995)(T = 25K at 40 AU). 

We try to constrain the disc mass by requiring the tidal 
arm and the disc orbiting star A to be optically thick. This 
requires an initial disc mass of 1.67 x 10~^Mq and a current 
disc mass of 1.05 x 10~^Mq. In order to estimate the latter, 
we have considered an outer disc radius of 57 AU, as esti¬ 
mated by Cabrit et al. (2006). This estimate is higher than 


that estimated by both Cabrit et al. (2006) and Osterloh Sz 
Beckwith (1995) {10~^Mq and 3 x lO^^M©, respectively), 
but lower than the ones estimated by Andrews & Williams 
(2005) (4 X 10~^ Mq). Note however that mass estimates 
from dust and CO isotopologues are always affected by very 
large uncertainties, and that our prescription for the tem¬ 
perature dependence is so simple that we do not expect to 
recover the observed mass with high precision. 

After adopting the temperature normalization and ini¬ 
tial disc mass as described above, the resultant line profiles 
show better agreement with the observations. In particular, 
both the A disc and the tidal arm are optically thick as ob¬ 
served, and their peak flux is in broad agreement with the 
observations by Cabrit et al. (2006), together with the mod¬ 
erate asymmetry of the double-horned line profile of the disc. 
The disc has in fact gained a small eccentricity, as apparent 
from the right bottom panel of Eig. 4. The ^^CO(J = 2 — 1) 
line profile of the gas around star B recovers the main fea¬ 
tures of the observed profile, besides the peak flux, which is 
lower than the observed one by a factor of 2. This may result 
from incorrect assumptions about the temperature profile 
around star B (temperature has a simple dependence on the 
distance from star A), but may also simply reflect that star 
B also possessed a disc prior to the encounter (in contrast 
to what we assume here, where the final disc is composed of 
material entirely captured from the disc of star A). 

Cabrit et al. (2006) obtained an indirect measure of the 
size of the disc around star A by using the position of the 
peaks of the typical double-horned profile of the emission 
lines (see Section 2.1). They observed a line profile with 
peaks at ~ 3.5 km/s from the centre, whereas in our syn¬ 
thetic observation we obtain the two peaks at ~ 2.5 km/s 
from the centre of the line. This seems to suggest that the 
disc in our simulations is larger than the one observed by a 
factor ~ 2. However, if we look at the hydrodynamical sim¬ 
ulation of our best model (see Eig. 9), or if we look at the 
synthetic continuum emission shown in Eig. 10, the disc size 
is in good agreement with the observations (~ 50 AU). After 
convolution we are obtaining contamination from low veloc¬ 
ity gas in the inner part of the tidal arm and the elliptical 
outer region of the disc. The difference in our synthetic line 
profiles implies that we are over-estimating emission from 
these outer regions. Detailed thermal modelling is however 
beyond the scope of this project; the agreement with the 
size derived from the thermal continuum reassures us that 
the dynamical effects of the encounter are well reproduced. 

5.2 Dust emission 

We used TORUS (c.f. section 3.2) to produce a synthetic ob¬ 
servation of the continuum emission at 1.3 mm of the same 
snapshot showed in Eig. 6. We consider standard assump¬ 
tions on the dust properties. Eirst of all we assume that 
spatial distribution of the dust follows exactly the gas dis¬ 
tribution. This is likely to be the case in our framework, 
since the dynamics of the gas itself is mainly ballistic during 
the encounter. The dust component is dynamically driven 
by the same gravitational potential. At the moment we ne¬ 
glect secondary effects, that might be due to hydro terms 
in the dynamics of the gas. We assume a gas to dust ra¬ 
tio equal to 100, and a grain size distribution ranging be¬ 
tween 50 A and 1 cm. We use a typical grain size distribu- 
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V (km/s) 


Figure 7. The integrated line profiles of the ^^CO(J = 1 — 0) emission line of the A disc (top left panel), tidal arm (top right panel) 
and molecular complex (bottom panel). The flux has been multiplied by a factor 4, in order to compare these synthetic profiles with 
the observed ones in Cabrit et al. (2006). In each panel the lines represent local thermal equilibrium (LTE) (blue), non-LTE with 
temperature profile estimated by Cabrit et al. (2006) (red) and non-LTE with temperature profile Infrared SED fitting by Osterloh & 
Beckwith (1995)(T = 25K at 40 AU, green line). Given that the three models are identical in all other respects, the LTE assumption 
indeed results in a lower integrated flux densities. 


tion, where dn/ds (x where n is the numerical density 
of dust particles, and s is the grain size. It is well known 
that in the interstellar medium w assumes a typical value 
w ~ 3.5 (Mathis, Rumpl & Nordsieck 1977). More recently, 
interferometric observations have suggested a shallower dis¬ 
tribution in the outer regions of protoplanet ary discs, where 
w 3 (e.g. Ricci et al. 2010a,b, 2012; Testi et al. 2014). 
Following these recent results, for our synthetic observation 
we use w = 3.3. We assume that the dust is comprised of 
amorphous silicate grains (Draine 2003a,b), obtaining an 
opacity = 0.I3Icm^/g at 1.3 mm (assuming a gas to 
dust ratio equal to 100). Fig. 10 shows the emission we ob¬ 
tain in mJy/beam, where the beam size is the one used by 
Cabrit et al. (2006) (0.87" x 0.56", see their Table I). Con¬ 
tour levels are submultiples and multiples of the noise level 
a = 0.47 mJy/beam in the same observation by Cabrit et al. 
(2006): (1/24, 1/6, I, 3, 48)cr. From disc A we obtain emis¬ 
sion levels that are compatible with the one observed (given 
the large uncertainties on the opacity), and a disc size that is 
in agreement with the one estimated by Cabrit et al. (2006) 


from the line profiles (~ 50 AU from the unconvolved emis¬ 
sion map). We underestimate the dust emission around star 
B, and this is again suggesting that our simple temperature 
treatment is underestimating the temperature of the mate¬ 
rial around the secondary star. By comparing the convolved 
emission map with the one shown by Cabrit et al. (2006) 
where the primary disc is located, it seems that we are over¬ 
estimating the total emission by a very small factor ~ 1.5 
(compare the size of the 3a contours), which is very small 
when compared to the uncertainties of the problem. Note 
that the normalization of the emission of our synthetic ob¬ 
servation does depend on our assumptions, mostly the dust 
opacity (e.g. in discs is often assumed a lower opacity at 
1.3 mm, as motivated by Beckwith et al. 1990), but also the 
uniform dust-to-gas ratio, the maximum grain size, and the 
grain size distribution index. With the assumptions we have 
made, the emission from the tidal tail is below the noise level 
(~ 0.3mJy/beam). 
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Figure 8. Line profiles of the synthetic observations of the best fit model for the disc around star A (top left panel), the tidal arm (top 
right panel) and molecular gas around star B (bottom panel). The red curve shows the ^^CO(J = 2 — 1) profile, and the blue curve the 
^^CO(J =1 — 0) one. The latter has been multiplied by a factor 4. The temperature profile is given by T{R) = 230 AT (A/AU)“°-^, and 
the mass of the disc is = 1.67 x 10 “^Mq. The synthetic profiles recover most of the features of the profiles derived by Cabrit et al. 
(2006), in particular flux intensities and kinematic signatures. The synthetic profile of the gas around star B is lower than the observed 
one by a factor of 2. This is due to the simple prescription we have used for the temperature radial profile, which fails to reproduce the 
emission properties of the gas close to star B, where the heating due to radiation from star B is certainly important. 


6 THE DIMMING EVENT 

Rodriguez et al. (2013) argued that the leading edge of the 
tidal arm may wrap around the A disc and obscure the line 
of sight between Earth and RW Aur A, thereby producing 
the observed dimming of RW Aur. In our simulations, there 
is indeed a stream of gas linking star A and B after the tidal 
encounter (however, it is not part of the primary tidal arm as 
Rodriguez et al. 2013, suspected, see the bottom left panel 
of Fig. 4). A portion of this stream of gas might provide the 
occulting body that caused the observed dimming event, as 
the integrated line-of-sight column density (~ 10“^gcm“^) 
exceeds the threshold column density (~ 7 x 10~^gcm“^) 
required to produce the dimming - assuming a flux decrease 
of the A star of 91% (Rodriguez et al. 2013) and a dust opac¬ 
ity at optical wavelengths k, = 3000cm^/g (Men’shchikov & 
Henning 1997). 

In order to verify this hypothesis, we have analysed the 
SPH simulation of the best-fitting model. We select all the 
foreground particles that are partially occulting star A at the 


present time by requiring h/D > 1, where h is the smooth¬ 
ing length associated with a single SPH single particle and 
D is the distance between the particle and star A in the 
plane of the sky. We then trace back their trajectories, in 
order to see whether the selected particles have a credible 
dynamical history. In Fig. 11 we plot all the SPH particles of 
the simulation with blue crosses, the two stars with yellow 
star symbols, and the occulting particles with filled dots. 
The orange dots represent the particles that show a credi¬ 
ble dynamical history. They start orbiting together within 
the circumstellar disc (see left-hand panels), and they then 
get incorporated in the well resolved bridge linking star A 
and B. At the present time they are located at the low den¬ 
sity periphery of the bridge (see right-hand panels). The 
occulting particles represented by the red dots are individ¬ 
ually scattered particles that originate from the very poorly 
resolved outer regions of the initial disc and are not asso¬ 
ciated with any coherent fluid structure. Since the orange 
particles are very few (only 4 for the particular choice of 
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Figure 10. Left panel: synthetic observation of our best-fit model in the continuum at 1.3 mm. Right panel: same emission map, convolved 
with the beam size used by Cabrit et al. (2006) (FWHM: 0.87" x 0.54"; PA = 20°). The colour code shows emission in mJy/beam. 
Contour levels are at (1/24, 1/6, 1, 3 and 48)f7, where cr = 0.47mJy/beam (see Table 1 by Cabrit et al. 2006). 



Figure 9. Circumstellar disc around star A as seen at t = 0 
(today) from the SPH simulation on the xy-plane. Column density 
is in code units. The outer radius of the disc is in agreement with 
values estimated by Cabrit et al. (2006) (40 — 57 AU). The disc 
has clearly gained an eccentricity after the encounter with star B. 

viewing angle), we cannot determine the structure of the low 
density margins of the bridge, and even less the size of the 
occulting body (0.09 — 0.27 AU, as estimated by Rodriguez 
et al. 2013). An accurate determination of the density struc¬ 
ture of the low density periphery of this bridge would re¬ 
quire an unfeasibly high resolution. However, the presence 
of particles in this region indeed indicates that possible ~ 
AU scale inhomogeneities at the edge of the stream linking 
A and B could cause the dimming event observed by Ro¬ 
driguez et al. (2013). The velocity in the plane of the sky of 


these particles ranges between 0.97 — 1.09 km/s and is thus 
consistent with the estimates from Rodriguez et al. (2013) 
(0.86 — 2.58 km/s). 

Finally, we have computed the evolution with time of 
the average column density lying between star A and the 
observer. In Fig. 12 we show that in the past the star has 
been completely occulted (at optical wavelengths) by the 
density structure linking stars A and B. This bridge has 
then evolved, until the line of sight only grazes it nowadays 
(t = 0). In the future, the average column density stays well 
below the value required to observe a dimming of star A as 
observed by Rodriguez et al. (2013). Moreover, for t ^ 0, 
the column density in the simulation starts to be dominated 
by particle noise. However, inhomogeneities in the bridge 
structure might entail new events in the next years. 


7 SUMMARY AND CONCLUSIONS 

We have used hydrodynamic models and synthetic observa¬ 
tions to constrain the observed morphology, kinematics and 
line fluxes of the binary system RW Aur which was first 
proposed as a candidate star-disc encounter by Cabrit et al. 
(2006). We draw the following conclusions from this work. 

With a suitable set of orbital parameters and viewing 
angle (See Table 1), we can reproduce almost all morpho¬ 
logical and dynamical features of the system: 

• RW Aur A has a circumstellar disc with an outer ra¬ 

dius of ~ 50AU which is within the observed range of 
40 — 57AU. The disc is optically thick in emission, 

as the ^^CO(J = 2 — 1)/ ^^CO(J =1 — 0) ratio is ~ 4. 

• RW Aur B is separated from star A by ~ 200 AU ( 
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Figure 11. 3D structure of RW Aur before (left panels, at —636 yr from today) and after (right panels, at 0 yr from today) the encounter 
for the best-fitting model. As in the other figures, the x-y plane represents the plane of the sky. The observer is located at 2 ; = + 00 . The 
origin of the axes is defined as the position of star A. Blue crosses represent all the SPH particles that have not been accreted on either 
of the two stars, which are shown in yellow. Filled dots indicate the particles that are partially occulting star A at present time. The red 
ones are scattered particles that do not show a credible dynamical history. The orange ones evolve in a coherent fluid structure. They 
are initially orbiting within the disc, are then stripped by the tidal interaction, and end up in the low density periphery of the stream 
connecting stars A and B. The dimming event observed by Rodriguez et al. (2013) can be explained by a clumpy inhomogeneity at the 
edge of such stream. Transverse velocities of the particles and column density are in broad agreement with the estimates by Rodriguez 
et al. (2013). 
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Figure 12. Average column density between the observer and 
star A as a function of time. The horizontal black line represents 
the column density that would cause a dimming of ~ 90%, as 
observed by Rodriguez et al. (2013). In the past the bridge con¬ 
necting stars A and B has been perfectly on the line of sight of 
star A, therefore completely occulting the star at optical wave¬ 
lengths. Nowadays such structure grazes the line of sight. Future 
dimming events cause by inhomogeneities in the density structure 
can not be ruled out by this model. 


~ 1.5'' for 140 pc). The projected AB separation is increas¬ 
ing at rate of ~ 0.002"/yr. The molecular complex on B is 
blueshifted with respect to star A with a broad distribution 
of radial velocity of 3 —6 km/s. The line profile from the star 
B is asymmetric and thus most likely due to material cap¬ 
tured during tidal encounter instead of a circumstellar disc. 
The molecular emission is marginally optically thick with a 
^2CO(J =2-1)/ 12C0(J = 1 - 0) > 4. 

• The expanding, spiral-shaped structure connected to 
RW Aur A disc is most likely a tidal arm produced in a 
fly-by of star B. The tidal arm extends about 600 AU ( ~ 4" 
for 140pc). The arm is entirely redshifted with a radial ve¬ 
locity ~ 3 km/s with respect to star A. This is faster than 
the escape velocity at corresponding distance, indicating the 
arm is indeed gravitational unbound. 

The synthetic molecular line data cubes obtained from a 
purely radial temperature estimated by Cabrit et al. (2006) 
have flux densities 3 — 5 times higher than observations. 
Nonetheless, by lowering the temperature to typical disc 
temperature for CTTS, we can reduce the flux densities to 
the observed values. 

The agreement between observations and simulations 
thus lends strong support to the hypothesis that the mor¬ 
phology of the RW Aur system is consistent with a tidal 
encounter scenario. Although the observational evidence for 
tidal encounters on Galactic scales is well established, this is 
the only system that is a good candidate for such encounters 
on the scale of protoplanetary discs (note that a tidal en¬ 
counter scenario has also been proposed as a possible model 
to explain the recent ALMA observations of AS 205 S, see 
Salyk et al. 2014). 

Finally we note that our simulations provide support 
for the hypothesis of Rodriguez et al. (2013) regarding the 
origin of observed dimming of star A in 2010/2011. We hnd 


that it is plausible that this was caused by inhomogeneities 
in the outskirts of a ‘bridge’ structure linking star A and 
B which may graze the line of sight at current epochs. Our 
simulations do not have the resolution to trace the structure 
of the low density periphery of the optically thick bridge, 
nor to model the ~ AU scale inhomogeneities that would 
be required to explain such an event. Nevertheless our sim¬ 
ulations predict a good agreement between the velocity in 
the plane of the sky of material in this region with that re¬ 
quired by the observed ingress time of the dimming event. 
This provides good support for the occultation hypothesis. 
It suggests that - depending on the hhing factor of such in¬ 
homogeneities - further dimming events should not be ruled 
out (as conhrmed by the observations by Antipin et al. 2015) 
and that these would be expected to share similar kinematic 
properties. 
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